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(54) System and method for the algebraic derivation of physiological signals 



(57) A technique for the algebraic derivation of a 
pulse oximetry signal uses the mathematical relations 
of detected light signals. First (6) and second (8) light 
sources transmit light through the patient's finger or re- 
flects light off the blood vessels in the patent's finger. A 
light detector detects light from each of the light sources 
and generates a measured intensity signal. The meas- 
ured intensity signal includes the true intensity of light 
transmitted from each of the light sources as well as 



noise introduced during the measurement process. A 
data sample from each of the light sources is digitized 
and subdivided into a plurality of data windows. The 
mathematical relationship (210) is apllied to each of the 
data windows (21 4) to determine a ratio value indicative 
of oxygen saturation in the patient. The ratio values (a, 
P) are statistically analyzed (224) to determine the best 
estimate for the correct ratio value (a). This ratio value 
is subsequently processed to determine the oxygen sat- 
uration within the patient. 
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Description 

Cross-Reference to Related Application 

s The present application is a continuation-in-part of United States Patent Application Serial Number 08/442,834, 

filed on May 17, 1995. 

Technical Field 

io The present invention relates generally to signal processing and, more particularly, to a system and method for 

processing physiological signals in the presence ol noise to derive the physiological signals. 



Background of the Invention 

15 The measurement of physiological signals is difficult because the underlying physiological processes generate 

very low level signals and interfering noise is inherent in the body and the interface between the body and sensors of 
the physiological processes. For example, the measurement of electrocardiogram (ECG) signals is based on the elec- 
trical activity generated by the electrical depolarization of the heart muscle. The signals are typically detected by surface 
electrodes mounted on the chest of the patient. The signals are initially weak at the signal source (/.a, the heart) and 

20 are even weaker at the surface of the chest. Furthermore, electrical interference from the activity of other muscles, 
noise caused by patient breathing, general movement, and the like cause additional interference with the ECG signal. 
External electrical interference, such as 60 Hertz (Hz) interference, also compounds the ECG measurement problem. 
Therefore, great care must be taken in the design and use of physiological processors to enhance the quality of the 
desired signal and reduce the effects of interfering signals. 

25 Another common physiological measurement that is made difficult by the presence of interfering noise is the meas- 

ure of oxygen saturation in the blood. This measurement is frequently performed with a pulse oximeter 1, illustrated 
in the functional block diagram of Figure 1 . A transmissive pulse oximetry sensor 2 is placed on a finger 4 of the patient. 
First and second light sources 6 and 8 are directed through the fleshy portion of the finger 4 and detected by one or 
more light detectors 10 on the opposite side of the finger. As is well known in the art, the light from light sources 6 and 

30 B are of different wavelengths that are differentially absorbed by oxygenated blood cells. The first light source 6 is 
typically designated as a Red light source having a wavelength in the red region of the spectrum. The second light 
source 8 is typically designated the IR source having a wavelength in the near infrared region of the spectrum. 

The pulse oximeter 1 determines the oxygen saturation based on a ratio of the light detected from the Red light 
source 6 and the IR light source 8, respectively. A ratio calculator 12 determines the ratio of detected light and uses 

3$ the value of the ratio as an address in a look-up table 14. The look-up table 14 contains data relating the ratio of 
detected light to the oxygen saturation in the blood. Atypical oxygen saturation curve 18 is illustrated in Figure 2 where 
the percentage of oxygen saturation is plotted against the ratio of detected light from the Red light source 6 and the 
IR light source 8 (see Figure 1). Pulse oximeters may also use reflective pulse oximetry sensors (not shown) in which 
the light sources and light detectors are positioned adjacent each other, and the light from the light sources is reflected 

40 back to the detector(s) by oxygenated blood cells in the finger 4. 

The measurement of blood oxygen saturation is important for physicians who are monitoring a patient dunng 
surgery and at other times. As with other physiological measurements, pulse oximetry measurement also is susceptible 
to interference from noise. As is known in the art, pulse oximetry is particularly susceptible to interference from stray 
light and from patient motion. Stray light detected by the light detector 10 can cause erroneous calculation of the ratio. 

45 Known techniques are employed to reduce the interference caused by stray light. The interference from patient motion 
is a much more difficult noise source and is the subject of intensive research. 

Therefore, it can be appreciated that there is a significant need for a system and method tor measurement of 
physiological signals that enhances the desired signal in the presence of interfering noise signals. This and other 
advantages provided by the present invention are described in the detailed description and accompanying figures. 



so 



Summary of the Invention 



The present invention is embodied in a system and method for the enhancenent of physiological signals in the 
presence ol noise. The system includes a plurality ol light sources directed toward the subject, each of which produces 
55 light of a different wavelength. A light detector is positioned to detect a plurality ol light signals after passage through 
or reflection Irom the subject and to generate signals indicative of the intensity of the plurality of detected light signals. 
Each of the detected light signals has a first portion arising from light transmitted from a corresponding one of the 
plurality of light sources and a second portion arising from a plurality of interference light sources. A window generator 



2 



EP 0 744 154 A1 



generates first and 6econd data windows derived from the plurality of detected signals. A storage location within the 
system contains a mathematical relationship of the first and second portions of the plurality of detected signals and a 
first ratio of the first portion of a first of the plurality of detected signals to the first portion of a second of the plurality of 
detected signals. A calculator uses the mathematical relationship to generate first and second functions using the first 
s and second data windows, respectively. The first and second functions provide a solution value indicative of the first 
ratio. 

The solution value may be a point of intersection between the first and second functions. Alternatively, the solution 
value is a statistical measure of the first ratio based on the solution value. 

In one embodiment, the first ratio is indicative of blood oxygen saturation in the subject. The system further includes 
10 a lookup table containing data relating the first ratio to the blood oxygen saturation level. The first data window is 
derived from a number of possible sources such as a derivative of the first of the plurality of detected channels. Alter- 
natively, the first data window may be the alternating current (AC) component of the first of the plurality of detected 
signals. 

In another embodiment, the system is capable of measuring carbon monoxide levels in the blood. In this embod- 

is iment, three light sources are directed toward the subject. The window generator generates first, second, and third 
data windows derived from the first, second, and third detected signals. The mathematical relationship also includes 
a ratio of the first portion of the first detected signal with the first portion of a third detected signal. The calculator uses 
the mathematical relationship to generate first and second functions whose solution value is indicative of a first ratio. 
The calculator also generates a third function and calculates a solution value of the second and third functions to 

20 determine a second ratio. The second ratio is related to blood carooxy^hemoglobin saturation. 

In yet another embodiment, the system generates a plurality of data windows, each providing a solution value 
indicative of the first ratio. A statistical analyzer analyzes the solution values and determines a final value for the first 
ratio based on the plurality of solution values. 

In one embodiment, the plurality of data windows are generated from a groupof signals comprising the first detected 

25 signal, a mathematical derivative of the first detected signal, the second detected signal, a mathematical derivative of 
the second detected signal, a combination of the first and second detected signals, and a mathematical derivative of 
the combination of the first and second detected signals. Alternatively, the system may calculate an average value for 
each of the group of signals and place data whose value is less than the average for a particular one of the group of 
signals in a first data window and placing data whose value is greater than or equal to the average for the particular 

30 one of the group of signals in the second data window. Alternatively, the second data window may be derived by placing 
data whose value is greater than or equal to the average for a second one of the group of signals in the second data 
window. 



Brief Description of the Drawings 

35 

Figure 1 is a functional block diagram of a prior art oximetry system. 

Figure 2 is a typical oxygen saturation curve employed by the system of Figure 1 to determine blood oxygen 
saturation. 

Figure 3 is a functional block diagram of a conventional adaptive signal processor. 
40 Figure 4 is a detailed functional block diagram of the system of Figure 1 . 

Figure 5 are waveforms that illustrate the liming control of light sources used by the system ol Figure 4. 

Figure 6 illustrates a waveform used in the calculation of a reference noise signal by the conventional adaptive 
signal processor of Figure 3. 

Figure 7 is a functional block diagram of the present invention used with the system of Figure 4. 
45 Figure 8 illustrates a first embodiment of the system of Figure 7. 

Figure 9 is a function block diagram of the window generation and analysis system of the present invention. 

Figure 1 0A is a graphical illustration of the equation used by the system of Figure 9 for a first data window. 

Figure 1 0B is a graphical illustration of the equation used by the system of Figure 9 for first and second data 
windows. 

so Figures 11 A-11C illustrate the use of multiple data windows by the system of Figure 9. 

Figures 12A and 12B are flowcharts illustrating the statistical analysis performed by the system of Figure 9. 
Figure 13A and 13B are flow charts illustrating one embodiment of the measurement process by the system of 
Figure 9. 

Figure 14 illustrates the use ol binned data used by the system of Figure 9. 
ss Figure 1 5 is a flowchart of an alternative technique used by the system of Figure 9. 
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Detailed Description of the Invention 

Measurement of physiological signals in the presence of interference is a difficult task, particularly if the interference 
is somewhat random rather than periodic. A number of different techniques can potentially be used to separate the 

s desired physiological signal from the interfering noise signal. For example, a filter can sometimes be used to remove 
the interfering noise signal. Notch filters, such as a 60 Hz notch filter, can be used to minimize interference from line 
noise Similarly, high frequency interference noise signals can be eliminated with a lowpass filter designed to pass the 
physiological signal of interest and to reject frequencies above the physiological signal bandwidth. However, some 
interference sources have the same or similar frequency content as the physiological signal of interest. For interference 

io of this type, different signal processing technologies must be employed. 

Adaptive signal processing is one well known technique for the separation of a desired signal from an interference 
signal. Adaptive signal processing is based on the assumption that the noise caused by the interference signal is 
uncorrected to the desired signal. A conventional adaptive signal processor, configured as a correlation canceller, is 
illustrated in the functional block diagram of Figure 3. An adaptive processor 20 has a signal input 22 and a noise 

15 reference input 24. The noise reference input 24 is fed to an adaptive filter 28. The adaptive filter 28 generates a filter 
output 30 that is subtracted from the signal input 22 in a conventional subtracter 34. The subtractor £4 generates an 
error signal 38 having a value designated herein as e that is fed back to the adaptive filter 28. The adaptive filter 28 is 
automatically adjusted so that the error signal 38 has a minimum correlation with the noise reference input 24. Thus, 
the adaptive filter 28 is adjusted so that the subtractor 34 cancels any correlated signal in the signal input 22. The error 

20 signal 38 is the system output and contains the portion of the input signal 22 that is uncorrected to the noise reference 
input 24 In a typical application of adaptive filtering, the signal input 22 consists of a combination of a pure input signal 
from a device, such as a sensor, and a noise signal from one or more sources. The noise reference input 24 should 
then be a signal that is related to and at least partially correlated with, the noise signal. The adaptive filter 28 is adjusted 
so that the error signal 38 is the pure input signal since the pure input signal has a minimum correlation with the noise 

25 reference signal applied to the noise reference input 24. 

Adaptive signal processing has been successfully applied to the measurement of physiological signals when the 
source of the interference signal is well characterized. For example, the physician may wish to listen to a fetal heartbeat 
whose acoustical signal strength is relatively small compared to the acoustical strength of the mother's heartbeat As 
discussed above, simple filtering will not work satisfactorily because the two heartbeats have similar frequency content. 

30 However adaptive signal processing can isolate the fetal heartbeat by using the much louder maternal heartbeat as 
the noise reference input 24 and the combination of fetal and maternal heartbeats as the signal input 22. Because the 
two heartbeats are uncorrected and the maternal heartbeat can be independently derived, the adaptive signal proc- 
essor 20 can easily isolate the fetal heartbeat Similarly, the adaptive signal processor 20 can remove 60 Hz interference 
by simply using the 60 Hz signal as the noise reference input 24. Thus, adaptive signal processing can effectively 

35 remove the undesirable interference signal provided that the interference signal can be independently derived. 

However, some physiological signals of interest do not have an independent interference source to use as the 
noise reference input 24. For example, pulse oximetry is susceptible to motion artifact, as described above. The motion 
alters the path that the light takes through the finger 4 (see Figure 1) and the characteristics of the interface between 
the finger 4 and the sensor 2. As the light from the Red light source 6 and the IR light source 8 pass through the fleshy 

40 portion of the finger 4, each is contaminated by a noise signal, primarily due to patient motion. The detected light is 
thus the combination of the true light transmitted through the finger 4 plus the interfering noise introduced in the meas- 
urement process. This may be illustrated by the following equations: 



45 



R=R* + N 0) 



r=r- + n (2) 

50 where R is the light intensity measured by the light detector 1 0 (see Figure 1 ), R* is the true intensity of light transmitted 
by the Red light source 6, and N is the noise source introduced by the measurement process while measuring the 
intensity of the Red light. Similarly, r in equation (2) is the light intensity measured by the light detector 1 0, r* is the true 
intensity of light transmitted by the IR light source B, and n is the noise source introduced by the measurement process 
while measuring the-intensily of the I R light. „ .. u . ^. . ■*. m 

55 The goal of the measurement process is to determine the ratio of the true intensity of Red light, R transmitted 

through the finger 4 to true intensity of IR light, r* transmitted through the finger. However, most pulse oximetry system 
determine the ratio of the measured signal (La. R/r) or some processed version of the measured intensities due to an 
inability to determine the true intensity. The ratio of intensities, whether it is the ratio of measured intensities, true 
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intensities, or some processed version of the measured intensities, is designated herein as r a . 

Some prior art pulse oximetry systems attempt to minimize the effects of motion artifact through conventional 
filtering or modulation of the intensity of the light sources 6 and B. However, these processing techniques are not 
particularly effective because the motion artifact is caused primarily by movement of venous blood in the tissues of the 

s finger 4 rather than from some external noise source such as stray light. Conventional filtering may remove some 
undesirable noise, but the frequency content of the motion artifact is similar to that of the desired signal. Modulation 
techniques may reduce interference from stray ambient light, but have little effect on motion artifact because the primary 
noise source (e.g., venous blood movement resulting from patient motion) originates in the measurement pathway. 
Thus, the ratio determined by many pulse oximetry systems is not accurate.- 

io it should be noted that the intensity of detected light varies with the patient's heart beat thus creating a time-varying 

pulsatile waveform. The pulsatile waveform contains an alternating current (AC) signal component and a direct current 
(DC) component. A more accurate determination of the ratio r a is given by the following equation: 



15 



25 



(Red/^/Redpc) (3) 

(' R AC /1R DC> 



where Red AC is the AC component of the intensity of the measured Red light, R, Redoc is the DC component of the 
intensity of the measured Red light, IR AC is the AC component of the intensity of the measured IR light, r, and IR^ is 
20 the DC component of the intensity of the measured IR light. In practice: the DC components tend to cancel each other 
out thus normalizing the resultant ratio of AC components. Thus, equations (I) and (2) above may be more accurately 
shown as: 



R(t) = R*(t) + N(t) (4) 

r(t) = r*(t) + n(t) (5) 

30 where R(t)=Red AC and r(t)=IR AC to reflect the time varying nature of the signals. 

A typical prior art transmissive pulse oximetry system 100 is illustrated in the functional block diagram of Figure 
4 where the sensor 2 contains the Red light source 6 and the IR light source 8, typically on the same side of the 
patient" s finger 4. The Red and IR light sources 6 and 8 are alternately activated by a timer 110. The activation timing 
of the first and second light sources 6 and 8 is illustrated in the waveform of Figure 5. The Red light source 6 is actrvated 

35 in the period T1 . Following the period T1 , the IR light source 8 is activated during the period 12. Following the period 
T2, neither the Red light source 6 or the IR light source B is activated during the period T3. The pulse oximeter uses 
the period T3 to detect stray ambient light and determine a baseline value to compensate for the stray ambient light. 
Compensation of stray light is well known by those of ordinary skill in the art and will not be discussed herein. The 
timer 110 repeats the pulsation of the Red light source 6 and the IR light source 8 in the manner described above. It 

40 should be noted that the intensity of the light from the Red light source 6 and the IR light source 8 is automatically 
adjusted by a closed-loop system to assure an acceptable detected signal level. This closed-loop gain control is well 
known in the art and need not be discussed herein. 

The detector 10 detects light transmitted through the fleshy portion of the finger 4. The signals generated by the 
light detector 10 are passed to a demultiplexer 112. The demultiplexor 112 is coupled to the timer 1 10 and is controlled 

45 by the timer 110 to generate an independent signal for the light detected from each of the light sources 6 and 8, 
respectively The time division multiplexing used by the system 100 is well understood and will not be discussed in 
detail herein. As discussed above, the timer 110 enables the Red light source 6 during the period T1 . During that same 
period T1 the timer also controls the demultiplexor 112 so that the detected signals from the Red light source 6 are 
routed to a data line 114. During the time period T2, the timer 110 enables the IR light source 8 and controls the 

so demultiplexor 112 so that the detected signals from the IR light source are routed to a data line 116. Each of the data 
lines 114 and 116 can be coupled to optional amplifiers 120. The amplified signals are coupled to the inputs of an 
analog to digital converter (ADC) 124 lhat digitizes the signal in a conventional manner. It should be noted that the 
amplifiers 1 20 may be integrally formed as part of the ADC 1 24. The ADC 1 24 may also include optional lowpass filters 
(not shown) to assure that the analog signals are bandlimited below Ihe Nyquist rate of the ADC. 

55 The demultiplexor 112 is shown as a separate component in Figure 4 for the sake of clarity. Those skilled in the 

art will recognize thai the demultiplexing function can also occur after the signal from the light detector 10 has been 
digitized. The present invention is intended to encompass all such conventional techniques for demultiplexing the 
signals from the light detector 10. 
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The ratio circuit 1 2 receives the digitized signals and uses the ratio of R(t)/r(t) to determine a location in the look- 
up table 1 4. Assuming that no motion artifact is present, the data entry in the look-up table 1 4 corresponds to the blood 
oxygen saturation. In reality, the ratio calculated by the ratio circuit 12 may be inaccurate because of the motion artifact. 
A technique has been developed to use the conventional adaptive signal processor of Figure 3 to eliminate the 
b motion artifact. A reference signal related to the motion artifact interference source is independently derived and applied 
as the noise reference input 24 to the adaptive signal processor 20 (Figure 3). The noise reference input 24 uses 
detected signals from the Red and IR light sources 6 and 8. These techniques are described in PCT Patent Publication 
No. W092115955, published on September 17, 1992. The system described in this publication generates a noise 
reference signal related to the interference noise and uses this noise reference in the correlation canceller version of 
io the adaptive signal processor 20 shown in Figure 3. The adaptive signal processor 20 uses the noise reference to 
cancel the noise in the measured signal thus resulting in a signal that is representative of the true signal {i.e., the 
measured signal minus the noise signal). 

The noise reference signal generated by the prior art pulse oximeter has the following form: 



is 



N(t) = R(t) - cor(t) (6) 



where N(t) is the time varying noise reference signal, R(t) is the time varying detected signal from the Red light source . 
6 {i.e., true intensity plus noise), r(t) is the lime varying signal from the detected signal from the IR light source 8 {i.e., 
20 true intensity plus noise) and co is a selected value of the ratio r B . Equation (6) has been empirically derived to model 
the noise source. 

As can be seen from Equation (6) above, the prior art pulse oximeter must determine a value for co in order to 
generate the noise reference signal N(t). As seen in Figure 2, the ratio of the light intensities and thus the value of co 
lies within a range from 0.5 to 3.0. The limitation in the range of values for co is imposed by the physiology That is, the 

25 oxygen saturation value lies between 100% and 0%, with the corresponding ratios lying between a value of 0.5 to 3.0, 
respectively. To compensate for variations in the sensitivity of the sensor 2, a range of ratio values from 0.3 to 3.0 is 
typically used. The prior art pulse oximeter takes advantage of the knowledge that the ratio must lie within the range 
from 0.3to 3.0 and scans the entire range of possible values for the ratio and inserts each of these values into equation 
(6) above. The noise reference signal for each possible value of the ratio r a is provided as the noise reference input 

30 24 (see Figure 3) to the adaptive signal processor 20. The adaptive signal processor 20 in turn generates the value e 
for each of the possible values of the ratio. A typical output of the value e versus the ratio r a is illustrated by a waveform 
48, shown in Figure 6. The best estimate of the value of w is given by a peak 50 or a peak 52 of the waveform 48. It 
is known that if the value of co corresponds to the peak 50. then N(t) in equation (6) equals C t n(t) where C 1 is a constant 
and n(t) is the noise source introduced by the measurement process while measuring the intensity of light from the IR 

35 source 8 (see Figure 5). If the value of co corresponds to the peak 52, it is known that N(t) in equation (6) equals C 2 r* 
(t) where C 2 is a constant and r*(t) is the true intensity of light transmitted by the IR light source 8. The value of co 
corresponding to the peak 50 is inserted into equation (6) above to generate a noise reference signal N(t) as the noise 
reference input 24 (see Figure 3) of the adaptive signal processor 20. The error signal 38 is the noise signal n(t) if the 
value of co corresponds to the peak 52. However, by selecting the value of co corresponding to the peak 50, the reference 

40 signal N(t) corresponds to the noise signal n(t). In this case, the correlation canceller adaptive signal processor 20 
cancels out the constant as well as correlated signals between the signal input 22 and the noise reference input 24 
such that the error signal 38 is the desired signal r*(t). The true output signals are provided to the ratio circuit 12 (see 
Figure 4) and processed in the manner previously described. 

The disadvantage of this prior art approach is that generating the value e for each of the possible values of :the 

45 ratio r a is a computationally difficult and time consuming approach to adaptive filtering in pulse oximetry. As those 
skilled 3 in the art can appreciate, real-time calculation of blood oxygen saturation is important to the physician. This 
real-time constraint can only be met with the prior art approach using expensive and powerful digital signal processor 
(DSP) hardware. 

An alternative approach to the measurement of blood oxygen saturation is described in U.S. Patent Application 
so Serial No. 08/442,834 and embodied in a system 1 80, shown in the functional block diagram of Figure 7, This approach 
provides a more efficient computational process that does not generate the noise reference required by the prior art 
approach. Rather, the desired signal {i.e., the true intensity) is directly generated and does not use correlation cancel- 
lation techniques in the adaptive signal processor. An analyzer 182 coupled to the ADC 124 (see Figure 4) receives 
digitized signals 1 84 representing the measured light intensity, R(t), from the Red light source 6, and digitized signals 
ss 1 86 representing the measured light intensity, r(t), from the IR light source 8. The analyzer 1 82 processes these signals 
using mathematical relationships between the measured signals and the true intensities, to generate a true intensity 
output 188 equal to the true intensity, R*(t), and a true intensity output 190 equal to the true intensity. r*(t). The math- 
ematical relationships are stored in a mathematical relationship storage area 190 for use by the analyzer 182. 
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The analyzer 1B2 generates the ratio r a ol true intensities (/.ft, R*(t)/r*(t)) in the process of generating the true 
intensity outputs 188 and 190. A ratio output 192 is coupled to the lookup table 14 to permit the determination of oxygen 
saturation in a conventional manner, "me output of the lookup table 14 is a value S p 0 2 corresponding to the blood 
oxygen saturation. The system 180 may also include an optional S p 0 2 peak detector 194 to generate signals indicative 

5 of the peak oxygen saturation. The true intensity outputs 188 and 190 are useful for monitoring the patient oximetry 
waveforms and for calculating continuous blood pressure measurements. Techniques for calculating blood pressure 
from pulse oximetry output waveforms are described in U.S. Patent No. 5,269,310. The advantage of the present 
invention is that the signal that is directly generated is the desired signal rather than the noise reference signal. Fur- 
thermore, the processing techniques of the present invention require far fewer computational steps, thus improving 

io the rate at which accurate data can be obtained. 

With respect to Figure 6, research has shown that the peak 50 corresponds to the ratio of the true intensities (/. 
ft, FT(t)/r*(t)), while the peak 52 corresponds to the ratio of noise intensities (/.ft. N(t)/n(t)). The following description 
provides details of the mathematical derivation of the reference signals representing the true intensities. For purposes 
of the following description, the ratio of the true intensities may be defined by the following equation: 



IS 



« = 53S . (7) 

r*(t) 



where a is the value of the ratio r a corresponding to the peak 50 (see Figure 6), R*(t) is the time varying true intensity 
20 of light transmitted from the Red light source 6 and r*(t) is the time varying true intensity of light transmitted from the 
IR light source B. The ratio of noise signals introduced by the measurement process is defined by the equation: 




(8) 



where p is the value of the ratio r a corresponding to the peak 52 (see Figure 6), N(t) is the noise introduced during the 
measurement of the light transmitted by the Red light source 6 and n(t) is the noise introduced during the measurement 
of the light transmitted by the IR light source 8. It is also known that the following constraint exists between a and p: 

0.3 < ct < P < 3.0 (9) 



because of the physiological nature of the signals. It is noted that the percentage of oxygen saturation is also a time- 
varying signal, changing by approximately 0.5% Over time. However, it is assumed that the blood oxygen saturation is 
35 constant over the short period required to perform the measurements. Thus, a and p can be considered ratio constants 
for purposes ol the present discussion. 

Given equations (4)-(5) and (7)-(8), it is possible to express the relationship between a and p using the following 
matrix equation: 



"1 


0 1 0" 




'R*(0" 
r*(0 




R(t) 


0 


I 0 1 






r(t) 


1 


-a0 0 


X 


N(t) 




0 


0 


0 i-p 




n(t) 




0 



(10) 



where it is assumed that a * p. As previously stated, it is known that the primary cause of noise in transmissive pulse 
oximetry measurements is motion artifact caused by the movement of venous blood in the finger 4. Thus, the value p 
in equation (8) is related to oxygen saturation in the venous blood. The assumption that a * p is based on the under- 
go standing that a is a measure of arterial blood oxygenation while p is related to venous blood oxygenation. As the body 
takes oxygen from the blood, blood oxygenation decreases as blood moves from the arterial portion of the circulation 
system to the venous portion of the circulation system. Thus, the measure of arterial oxygenation, measured by a, is 
not the same as p, which is related to venous oxygenation. 

The significance of equation ( 1 0) is that all signal components can be explicitly calculated as a function of the input 
55 signals and the ratio constants a and p. The true signal components, R*(t) and r*(t) can also be explicitly derived using 
equation (10) above. The true signal components, R*(t) and r*(t), can be expressed in terms of the measured signals, 
R(t) and r(t), by the following equations, which are derived from equation (10): 
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15 



aR(t)-« P r(t) 
v ' a-p 

r(t) = 0CLiMQ (12) 

Similarly, the noise signals, N(t) and n(t), can be expressed in terms of the measured signals. R(t) and r(t). by the 
following equations, which are also derived from equation (10): 

N(t)= «Pr(t) p PR(t) (13) 

n{l)= «£flLRO (14) 



It will be noted that the above equations (11)-(14) provide the true signal components, R*(t) and r # (t), and the noise 
components, N(t) and n(t), as a function of the measured signals, R(t) and r(t), available from the sensor 2 (see Figure 
4) and the ratio constants a and p. The values of the ratio constants a and p are not known and must be determined 
by the system 1 80. 

20 The measurement system described in PCT Patent Publication Number W092115955. published on September 

1 7, 1 992, requires a determination of a peak value for the ratio r a . This system steps through each and every possible 
value of the ratio r a over the physiological range for the ratio r a to determine the peak value. Only then can the peak 
value be substituted into equation (6) to generate the noise reference signal. This is a tedious and time consuming 
process to determine the blood oxygen saturation. The correlation cancellation techniques described above with re- 

25 S pect to PCT Patent Publication Number W092115955 require a powerful and expensive digital signal processor to 
determine the blood oxygen saturation in real-time; 

Various embodiments of the analyzer 182 are described in U.S. Patent Application Serial No. 08/442,834. The 
analyzer 182 does not require a noise reference signal generated by the measured signals as does the prior art oxi- 
meter. Rather, the analyzer 182 directly derives a true intensity output 188 corresponding to the true intensity R* of 

30 light transmitted through the finger 4 from the Red light source 6 (see Figure 4) and a true intensity output 1 32 corre- 
sponding to the true intensity r* of light transmitted through the finger from the IR light source 8. The system 180 uses 
the ratio of R*(t)/r*(t) {Lo., a) and the oxygen saturation curve 1 8 (see Figure 2) to determine the blood oxygen saturation 
in a conventional manner. 

A first embodiment of the analyzer 1 82, shown in the functional block diagram of Figure 8, uses an adaptive signal 
35 processor 210. Although similar to the adaptive signal processor 20 of Figure 3, the adaptive signal processor 210 
does not use correlation cancellation techniques with a noise reference signal. Rather, the adaptive signal processor 
has an adaptive filter 21 2 with a filter output 214 that directly generates the desired output signal R*(t) if the appropriate 
signal is selected for a reference input 216 to the adaptive filter. 

A subtractor 220 has a positive subtractor input 224 and a negative subtractor input 228. The measured signal R 
40 (t), which is the combination of the true signal. R*(t), and the noise signal, N(t), is coupled to the positive subtractor 
input 224, while the filter output 214 is coupled to a negative subtractor input 228. The subtractor 220 generates an 
error signal 232 that is fed back to the adaptive filter 212 in a well known manner. The adaptive signal processor 210 
uses an iterative process to adjust the adaptive filter 212 to minimize the error signal 232. Minimization techniques, 
such as least mean squares (LMS) or least squares lattice (LSL), are used to adjust the adaptive filter 212. These 
45 techniques are well known in the art of adaptive signal processing and need not be discussed herein. 

The reference input 21 6 is provided with a signal R'(t) derived from equation (1 1 ) to estimate the true intensity R* 
(t). The signal R'(t) is simply the signal of equation (11) for selected values of the ratio r a over the range from 0,3 to 
3.0 to determine values for the ratio r a corresponding to the peaks 50 and 52, respectively The analyzer 1 82 does not 
scan the entire range from 0.3 to 3.0 as does the prior art pulse oximeter. In contrast, only selected values for the ratio 
so r a between 0.3 and 3.0 are used to determine the correct values of the ratio constants a and p thus resulting in a more 
computationally efficient approach to pulse oximetry. Furthermore, the prior art reference signal of equation (6) must 
be used as a reference signal in the correlation cancellation adaptive signal processor 20 of Figure 3. so that the error 
signal 38 is the desired signal. In contrast, the analyzer 182 of the present invention directly generates the desired 
signals using the mathematical relations of equation 10. When the correct values for the ratio constants a and p have 
been determined, the function R'(t) = R'(t). Again, it should be noted that the signal generated by the analyzer 182 is 
mathematically derived and equals the desired true intensity if the correct values are selected for a and p. This approach 
is markedly different from the prior art approach to adaptive signal processing because no noise reference signal is 
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generated and no noise canceller is used by the adaptive signal processor 210. The true signal is determined directly 
from the given conditions and the mathematically derived relationships shown in the equations above. The adaptive 
filter 212 can be designed in a well known manner to improve the accuracy and correctness of the true signal. The 
procedure for the selection of the proper values for the ratio constants a and p is discussed below. 

s (t should be noted that the above discussion relates to the measurement of the true intensity of light transmitted 

from the Red light source 6. However, those skilled in the art can readily recognize that the same principles apply to 
the measurement of the true intensity of light transmitted from the IR light source 8. The true intensity signal r*(t) can 
be directly derived from the true intensity signal R*(t) using the relationship of equation (7). Thus, both true intensity 
signals R*(t) and r*(t) can be directly derived once the correct values have been determined for the ratio constants a 

io and p. 

As stated above, the signal R'(t) provided to the reference input 216 is equation (11) for selected values of the 
ratio r a . The system 180 determines values for a and p so that R*(t)=R*(t), to assure that the filter output 214 will 
represent the true signal intensity R*(t). As previously illustrated by equation (9) above, the value of the ratio constants 
a and p lie between 0.3 and 3.0. The system 180 uses a peak detector 250 to derive the values of the ratio constants 

« a and p without scanning the entire range. The peak detector 250 detects a peak in the error signal 232 over the 
physiological range of values for the ratio r a . The physiological range for the ratio r a is divided in half and the peak 
detector 250 determines whether a peak lies in the first half of the subdivided range or the second half. If the peak 
detector 250 detects a peak in the first half of the subdivided range for the ratio r a , it discards the second half. Conversely, 
if the peak is detected in the second half of the physiological range, the first half is discarded. The peak detector 250 

20 repeats this process (i.e., divide and peak detect) until a peak is detected with sufficient accuracy. This approach 
provides a great computational advantage over the prior art since lar fewer calculations are performed to detect the 
peak value. 

The analyzer 182 (see Figure 7) produces the ratio output 192, and the value for oxygen saturation S p 0 2 may be 
determined in a conventional manner. The optional peak detector 194 may be used to determine peak S p 0 2 levels. 

2S Thus, the analyzer 812 (see Figure 7) directly produces reference signals equal to the true intensities. In practice, 
these true intensity signals are derived from the filter output 214. This direct calculation of the true intensities is per- 
formed without having to generate a noise reference signal as is done in the prior art, and without having to use digital 
signal processing correlation cancellation techniques that require a significant number of computational steps. Fur- 
thermore, the analyzer 1 82 requires significantly fewer calculations to determine accurate values for the ratio constants 

30 a and p. 

* Other embodiments of the analyzer 182 are also illustrated in U.S. Patent Application Serial No. 08/442,834. How- 
ever, these embodiments need not be described herein. The pulse oximetry system described in U.S. Patent Application 
Serial No. 08/442,834 is computationally efficient because it rapidly determines the peak value for the ratio r a without 
having to scan the entire physiological range of values as do systems of the prior art. However, this process still requires 

35 a powerful computational device such as a digital signal processor to generate the adaptive filter. Those skilled in the 
art will readily recognize that the cost of a typical digital signal processor can make a product undesirably expensive. 

The techniques described in the present invention generate the desired results without having to determine the 
peak values for the ratio r a and, furthermore, without requiring an expensive digital signal processor Indeed, the prin- 
ciples of the present invention can be implemented on a conventional digital computer. 

40 The present invention is embodied in a system 300, shown in the functional block diagram of Figure 9. The digitized 

signals 184 (see Figure 7) representing the measured light intensity R(t), from the Red light source 6, are stored in a 
first data storage area 204, The digitized signals 186 representing the measured light intensity, r(t), from the IR light 
source 8 are stored in a second data storage area 206. The first and second data storage areas 204 and 206 store 
the digitized signals 184 and 186, respectively, for a predetermined measurement period. The predetermined period 

45 is selected such that the oxygen saturation value can be considered a constant. In one embodimenVthe first and 
second data storage areas 204 and 206 store approximately five seconds of the digitized signals 184 and 186, respec- 
tively. 

A mathematical processor 210 receives the digitized signals from the first and second data storage areas 204 and 
206. The mathematical processor 210 calculates AC values and derivatives of the digitized signals 184 and 186. As 

so will be described in detail below, the mathematical processor 21 0 may also calculate squared values of digitized signals 
1 84 and 186 and calculates averages of the squared signals, the AC signals, and the derivatives. A window generator 
214 separates the digitized data into a plurality of pairs of windows based on the characteristics of the data and certain 
measurement criteria. The characteristics of the data in the measurement criteria used by the window generator 214 
will be described in detail below. 

55 A calculator 21 6 calculates a plurality of (a, P) pairs for each of the plurality of pairs of windows determined by the 

window generator 214. Thus, the output of the calculator 21 6 is a set of a and p pairs. The alpha values are stored in 
an a value storage area 21 8, and the p values are stored in p value storage area 220. A statistical analyzer 224 receives 
the a and p value data from the a value storage area 2 1 8 and the p value storage area 220, respectively: The statistical 
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analyzer 224 performs statistical analysis on the a and p value data to determine the most likely value for a. The output 
of the statistical analyzer 224 is the correct value a for the digitized signals 184 and 186. The oxygen saturation and 
peak oxygen saturation may be determined using the lookup table 14 (see Figure 4) and the peak detector 136 in a 
conventional manner. The analysis steps performed by the mathematical processor 210, the window generator 214, 
the calculator 216. and the statistical analyzer 224 do not require the sophisticated and powerful data processing 
capabilities of a digital signal processor, "merefore, the system 200 may be implemented on a conventional digital 
computer and is capable of producing real time results at a significantly lower cost than systems employing digital 
signal processor technology. 

It can be shown that the ratio constants a and p are interrelated. II one assumes that the true signal and the noise 
signal are uncorrected, the signals r*(t) and n(t) are said to be 'orthogonal." This may be defined mathematically by 
the following equation: 



JY*(t)n(t) = 0 



which may also be expressed as: 



(15) 



(ar (t)-R(t))(R(t)-pr(t)) _ 0 

(«-P) 2 (16) 



by conventional mathematical manipulation of equations (10)-(15). 

Using the previous assumption that a * p. the following equations relate the ratio constants a and (3: 

jR 2 (t)-p jR(t)r(t) 

a= L_ (17) 

jR(t)r(t)-p Jr-(t) 



jR 2 (t)-a jR(t)r(t) 

P= J - ? l —f~2 (18) 

jR(t)r(t)-a Jr 2 (t) 



which solve equation (16) for a and p, respectively. 

As seen in equations (17) and (18), the ratio constants a and p are symmetric and thus only one independent 
variable, either a or p, need be determined. The following description provides an example of the determination of the 
values of the ratio constants a and p. t 

As discussed above, the ratio constant p is related to oxygen saturation in the venous system. While a curve similar 
to that of Figure 2 has not been developed to indicate the oxygen saturation for the venous system, it is known that 
the ratio constant p can provide some measure of oxygen saturation in the venous system. For purposes of the present 
invention, it is assumed that oxygen consumption in the tissue is constant over the short duration of the measurement 
process. 

There are an infinite number of (a.p) pairs that provide valid solutions to equations (17) and (18). The equations 
are illustrated graphically in Figure 10A where the digitized data 184 and 186 provide a set of valid (a,p) pairs. An a 
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curve 260 provides a graphical representation of all valid values (or a and p that solve equation (17). Similarly, a p 
curve 264 provides a graphical representation of all valid values for a and p that solve equation (18). The symmetry 
of the a curve 260 and the p curve 264 graphically illustrate the symmetry between a and p previously discussed with 
respect to equations (17) and (18). 
5 if the same equations (17) and (18) are solved for a second data window, the result would be a second set of valid 

(a,p) pairs. The second data window is generated by the window generator 214 (see Figure 9) from the same digitized 
signals 184 and 186 used to generate the first data window. The calculator 214 generates the second set of valid (a, 
P) pairs, which are graphically illustrated in Figure 10B. The techniques used to generate the data windows will be 
discussed in detail below. 

io A second a curve 260' provides a graphical representation of all valid values for a and p that solve equation (17) 

for the second data window. Similarly, a second p curve 264' provides a graphical representation of all valid values for 
a and p that solve equation (18) for the second data window. As discussed above, the system 200 assumes that the 
oxygen saturation remains constant over the short time required to collect the data for analysis. Thus, the data from 
the first data window and the data from the second data window should result in the same values for a and p. In this 

is ideal setting, there is only one correct value for a and one correct value for p. The correct values are indicated by the 
intersection of the curves shown in Figure 10B. Specifically, an intersection 262 is the intersection of the first and 
second a curves 260 and 260' and is the only value for aand p that satisfies equation (17) for both the first and second 
data windows. Similarly, an intersection 266 is the intersection of the first and second p curves 264 and 264' and is the 
only value for a and p that satisfies equation (18) for both the first and second data windows. The symmetry of equations 

20 (17) and (18) is also illustrated by the intersections 262 and 266, which have identical coordinates, with the ordinate 
and abscissa values interchanged. 

While measurement technique used by the system 200 is illustrated graphically in Figure 10B, the system uses 
an analytical approach to solving for the correct values of a and p. Equation (17) may be rewritten in the form illustrated 
below: 



25 



30 



35 



40 



45 



o,=F,(P) 

a 2 = F 2 (P) (19) 

so as to express a as a function of p for the first and second data windows, respectively. The specific form of equation 
( 1 7) is rewritten as follows : 

jR?(0-p-jRi(t)Tj.(t) 
aj = Fi(P) = -L, ! j- 2 (20) 



where i = 1 to 2 for the two data windows, respectively. Since, for a short period of time a and p are constant, oh = otg. 
Thus, it is possible to set F 1 (P) = F 2 (p) and solve for a and p. The solution yields a pair of quadratic equations having 
the following form: 



so 
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Jr, 2 (t)fRl(t)-jR?(t)fr 2 2 (t) 



t t 



frf(OjR!(t)-Jr, a (t)jRi(t) 
t t 



\2 
) 



-4x Jr, 2 (t)jR,(t)r 2 (t)- fR,(t)r.(t)Jr 2 2 (0 
» t t J 

\ 

jR,(0r,(t)fR|(t)-jR?(t)jR2(t)r2(t) 



t t 



1/2 



-i- 1 



2x 



f r, 2 (t) JR 2 (t)r 2 (t)- JR , (t)r, (l) Jr| (t) 

W t t I 



(21) 



This process results in a pair of values for (o,P) so long as different data windows are used for F^P) and F 2 (P). 

The selection of data windows by the window generator 214 (see Figure 9) is illustrated in Figures 11 A-11C. In 
Figure 1 1 A, a data sample 300 is provided to the system 200 for each of the digitized signals 1 84 and 1 86. As previously 
noted, the precise length of the data sample 300 is not critical. However, the data sample 300 must be sufficiently short 
in length to validate the assumption that oxygen saturation is constant. In one embodiment, the data sample 300 may 
be five seconds in length. However, the principles ol the present invention are not limited by the precise length of the 
data sample 300. In Figure 11B, a first data window 302 and a second data window 304 are derived Irom the data 
sample 300. In the case of the two data windows 302 and 304, the system 200 will provide two (a, P) pairs. As illustrated 
in Figure 10B, the correct a and p, values may be ascertained by the intersections 262, and 266, respectively. As 
illustrated in Figure 11C, first, second, and third data windows 306, 308, and 310 result in three (a. P) pairs. The first 
(a, p) pair results from solving equation (21) using the first and second data windows 306 and 308. The second (a, p) 
pair results from solving equation (21) for the second and third data windows 308 and 310 while the third (a, P) pair 
results from solving equation (21 ) for the first and third data windows 306 and 310. Thus, any combination of two data 
windows from the data sample 300 will provide an (a, P) pair. It should be noted that the data windows will frequently 
overlap in time, such as illustrated in Figure 11 C. The data windows may be generated in a variety of different manners. 
For example, the data window 308 may be the derivative of the data sample 300, while the data window 310 may be 
the AC component of the data sample. 

The window generator 214 (see Figure 9) may generate data windows in a number of different manners. In one 
embodiment, only two data windows are generated. For example, the data windows may be diflerenl time segments 
of the data sample, such as illustrated by the first and second data windows 302 and 304 in Figure 11 B. Alternatively, 
the first data window 302 is generated using the AC component ol the digitized signals 184 and 1 86 for the entire data 
sample 300, while the second data window 304 is generated by taking the derivative of the digitized signals for the 
entire data sample. It should be noted that whenever a data window is generated for the digitized signal 184, the 
corresponding data points from the digitized signal 186 form part of the same data window. For example, consider a 
data sample in which the digitized signals 184 and 186 are each 500 data points long. Whenever a data point from 
one of the digitized signals 184 and 186 are placed in the first data window, the corresponding data point from the 
other of the digitized signals is also placed in the same data window. In the example above, if data point n, where n = 
1 to 500 in the above example, of the digitized signal 184 is placed in the first data window, then data point n of the 
digitized signal 186 is also placed in the first data window. Conversely, if a particular data point, such as data point n+1 
of the digitized signal 184 is placed in the second data window, then data point n+1 ot the digitized signal 186 is also 
placed in the second window. Thus, whenever first and second data windows are generated by the window generator 
214, it contains the corresponding data points from both digitized signals 184 and 186. In this manner, equation (20) 
uses corresponding data points from the digitized signals 184 and 186 in each of the respective data windows. 
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To determine the derivative of the digitized signals 184 and 186, the mathematical processor 210 (see Figure 9) 
performs a point-to-point derivative in which the elope from data point n to data point n+1 is calculated tor each of the 
data points in the first data storage area 204. A similar calculation is performed on the digitized signal 1 86 for the data 
points stored in the second data storage area 206. The AC component of the digitized signals 1 84 and 1 86 are deter- 

5 mined by calculating an average value for all data points in the first data storage area 204. This average value is 
subtracted from each of the data points in the first data storage area. A similar calculation is performed on the digitized 
signal 186 and the data points stored in the second data storage are 206. This calculation effectively removes the DC 
component of the digitized signals 184 and 186. 

Alternatively, the window generator 21 4 (see Figure 9) generates data windows based on both statistical knowledge 

to and clinical knowledge about the nature of the physiological signal. The mathematical process 210 (see Figure 9) 
calculates six average values from the initial data sample 300. These six average values are the square of the AC 
component (Red AC 2 ) of the measured signal R(r) from the Red light source 6 (see Figure 4), the square of the derivative 
of the measured signal from the Red light source (dR(t)/dt 2 ), the square of the AC component (IRac 2 ) of measured 
signal from the IR light source 8, the square of the derivative of the measured signal from the IR light source (dr(t)/ 

is dt 2 ) the product of the AC components of each of the measured signals from the Red light source and the IR light 
source (Red AC * I R AC ) and the derivative of the product of the measured signals ((dR(t)/dt) * (dr(t)/dt)). As those skilled 
in the art will readily recognize, calculating the average value of the AC component itself would result in a value of 
zero. For this reason, the average values calculated by the mathematical processor 21 0 are typically squared values. 
The window generator 214 selects data points for each data window based on the size of a particular data point 

20 compared to the average value tor that particular parameter. For example, one pair of data windows may be derived 
based on the measured signal from the Red light source 6. In this example, all'data points for which Red AC is less than 
or equal to the average AC component will be in a first window, while all data points for which the Red AC is greater 
than or equal to the average value will be in a second data window Similar comparisons are performed for each of the 
remaining five averages described above. Thus, for each of the six averages calculated by the mathematical processor 

25 21 0 the window generator 21 4 divides the data sample 300 into a pair of windows with data points that are less than 
or equal to the average for a particular parameter being placed in one data window for the parameter, while the data 
points that are above average for that parameter are placed into a second data window 

In addition to the six pairs of data windows derived from the six averages described above, the window generator 
214 generates two additional data windows based on the relative size of the AC components of the measured signals 

so (Red AC and IR AC ). If a particular data point for the Red AC signal is less than or equal to the corresponding data point 
tor the IR AC signal, that data point is assigned to a first data window. If the particular data point of the Red AC signal is 
greater than or equal to the corresponding data point of the IR AC signal, that data point is placed in the second data 
window Similarly, an additional pair of data windows is generated by the window generator 214 based on the relative 
sizes of the derivatives of the measured signals (dR(t)/dt) and dr(t)/dt). Thus, the window generator 21 4 can generate 

35 eight pairs of data windows from the data sample 300. The optimum window pairs would have all noise in one window 
and all signal in the corresponding window. However, the data points within a particular window are generally a mixture 
of both noise and signal. However, based on the characteristics of the data, and the placement of data in one window 
or another based on its value relative to an average value for a particular parameter, it can be expected that, for each 
pair of data windows, one data window will contain more signal, while the other corresponding data window will contain 

40 more noise. . . 4U ... 

As previously discussed, if a particular data point derived from the digitized signal 184 is placed in the first data 
window, the corresponding data point derived from the digitized signal 186 is also assigned to the first data window. 
In contrast, H the particular data point derived from the digitized signal 184 is- placed in the second data window, the 
data point the window generator would 214 also assign the corresponding data point derived from the digitized signal 

45 1 86 to the second data window regardless of its actual value. In this manner, a pair of data windows may be generated 
based on the average value of a particular signal derived from the digitized signal 184. 

Alternatively, the window generator 214 may analyze the data points derived from the digital signal 186 and classify 
those data points as belonging in the first or second data window based on their value relative to the average value of 
the particular parameter derived from the digitized signal 186. If the window generator 214 performs this analysis on 

so the data points derived from the digitized signal 186. the corresponding data points derived from the digitized signal 
184 are assigned to the first or second data window based only on the analysis of the data points derived from the 
digitized signal 1 84 and not based on the value of the data points derived from the digitized signal. 

The calculator 216 (see Figure 9) may use any number of combinations of the windows generated by the window 
generator 214. For example, the calculator 216 may use the AC components, R AC and r AC , and the derivative values 

ss dR(t)/dl, and dr(t)/dt, to generate six diflerent window pairs, as shown in Table 1 below. 
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As those skilled in the art can appreciate, a large number of possible data window combinations can be generated by 
the window generator 214 using the average values discussed above. The present invention is not limited by the 
specific combination of data windows used by the calculator 216. As discussed above, the average value calculated 
by the mathematical processor is used by the window generator 214 as a threshold value. Any data points in the first 
data storage area that are less than or equal to the threshold value are placed in the first data window, while data 
points that exceed the threshold value are placed in the second data window. As discussed above, the corresponding 
data points in the second data storage area 206 are placed in the same data window as the corresponding data point 
in the first data storage area 204. By selecting an average value, and classifying data points based on this average 
value, it can be expected that for each pair of data windows, one data window will contain more noise while the other 
will contain more signal. 

As discussed above, equation (21) results in a pair of valid (<x,P) pairs under predetermined conditions. It is not 
initially known which of the values corresponds to a and which value corresponds to p. However, the value for a can 
subsequently be determined imposing the following clinical condition: 



0-3 * Roc/roc < a < p < 3 0 * Roc/roc 



(22) 
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It should be noted that equation (22) is identical to equation (9) above except for the introduction of the DC component 
ratio in equation (22). The DC components merely indicate that the values for a and P have not been normalized in 
the manner discussed with respect to equations (3)-(5). As those skilled in the art can appreciate, normalization of the 
measured signals, R(t) and r(t), prior to analysis by the system 200 could result in the loss of useful data prior to 
calculating the average values discussed above. 

The calculator 216 (see Figure 9) solves equation (21) for a and p. The values of a and p determined by the 
calculator 216 typically satisfy equation (22). However, in some circumstances, the values of a and p do not meet the 
requirements of Equation (22). For example, in certain circumstances, the digitized signals 1B4 and 186 may contain 
significant amounts of noise. 

The statistical analyzer 224 (see Figure 9) uses statistical techniques to determine an accurate value for a. The 
above example illustrated the operation of the system 200 using only two data windows generated by the window 
generator 214. In one embodiment, the system 200 calculates a new (a, P) pair approximately once per second. This 
new (a, p) pair is analyzed along with the previous 11 (a. p) pairs to drive a statistical measure of a. This process is 
illustrated in the flowchart of Figures 12Aand 12B. 

At the start 312. the system 200 has a new sample of the digitized data 184 and 186 (see Figure 9). In step 314, 
the statistical analyzer 224 calculates the average values for a and p from 1 2 (a, p) pairs including the present (a. p) 
pair and the previous 11 (a, p) pairs. In addition, the statistical analyzer 224 calculates the mean deviation for the 12 
(a, P) pairs. In decision 316, the statistical analyzer determines whether the deviation of a is less than 0.05. If the 
deviation of a is less than 0.05, the system accepts the new value for aand ends the process in step 352. If the deviation 
of a is not less than 0.05, the result of decision 316 is NO. 

In that event, the statistical analyzer 224 performs a measurement on each of the 12 (a, P) pairs. For each of the 
12 (a, P) pairs, the statistical analyzer tests, in decision 318, whether p - a is greater than 2.0. If p - a is greater than 
2.0. the result of decision 318 is YES. In that event, one value of the particular (a. p) pair is considered invalid. In step 
322, the statistical analyzer keeps the value that is closer to either the a or the p average. This process is repeated 
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lor each of the 12 (a, P) pairs. 

Following the analysis of each of the 12 (a ( P) pairs, the statistical analyzer 224 enters a loop in which data is 
processed until the deviation of a is less than a predetermined value. In step 326, the statistical analyzer 224 recalcu- 
lates average values for a, P, mean deviation of a, and mean deviation of p. In decision 330, the statistical analyzer 

s 224 determines whether the deviation of u is less than 0.05. If the deviation of a is less than 0.05, the result of decision 
330 is YES. In that case, the present value for a is accepted, and the system ends the process in step 352. If the 
deviation of a is not less than 0.05, the result of decision 330 is NO. In that event, the statistical analyzer 224 performs 
additional analysis tor each of the 1 2 (a, P) pairs. In decision 332, the statistical analyzer 224 tests to determine whether 
the value for p is closer to the average value for a than it is to the average value for p. If the value lor p is closer to the 

w average a value than it is to the average p value, the result of decision 332 is YES. In that event, in step 334, the 
statistical analyzer 224 discards the value for a and sets a equal to the value for p. Following step 334, the statistical 
analyzer 224 moves to decision 336. If the p value is not closer to the average a value than it is to the average p value, 
the result of decision 332 is NO. In that event, the statistical analyzer moves to decision 336 where it determines 
whether the value for a is closer to the average p value than it is to the average a value. If the a value is closer to the 

is average p value than it is to the average a value, the result of decision 336 is YES. In that event, in step 340, the 
statistical analyzer 224 discards the value for p and sets p equal to the value of a. If the value of a is pot closer to the 
average p value than it is to the average a value, the result of decision 336 is NO. In that event, or upon completion 
of step 340, the statistical analyzer 224 repeats this process lor each of the 1 2 (a, p) pairs. 

Alternatively, the statistical analyzer 224 determines the confidence factor, which is a measure of the confidence 

20 in the calculation ol a. As previously described, there are certain circumstances in which the calculated value for a is 
discarded and the value for p is used in its place {i.e., p is relabeled as a). To determine the confidence factor, the 
statistical analyzer 224 tracks the previous 1 2 a values and determines the average value, and the average deviation. 
If the calculated value for a is far from the average deviation of the previous 1 2 measurements, and the current value 
for p is approximately equal to the average value, the statistical analyzer 224 discards the value for a and uses the 

25 value for p in its place. The statistical analyzer 224 then determines the confidence value based on how many calcu lated 
values were discarded in the previous 12 measurements. If no a values were discarded in the previous 12 measure- 
ments, the confidence factor is 1 2/1 2 = 1 00%. However, if one a value was discarded in the previous 1 2 measurements, 
the confidence factor is 11/12 = 92%. In this manner, the system 200 provides a statistical measure of the confidence 
in the current measurement for a. 

30 Another potential drawback to the measurement technique used by the system 200 is possible system instability. 

With reference to Figure 10A, it is known that the first a curve 260 and the first p curve 264 asymptotically approach 
an intersection 270 between the first a curve and the first p curve. If there is any error in the value one parameter in 
the area of the intersection 270, the error is magnified for the second parameter. For example, consider two values for 
a, designated as and Og, respectively, on the first p curve 264. Although the values for and differ by less than 

35 0.1 . this difference results in a potential error of more than 0.3 in the value of p. This type of error magnification can 
lead to potential system instability. However, the use of multiple windows avoids the system instability by assuring that 
accurate values are selected for a and p. 

Following this process for each of the 12 (a, P) pairs, the statistical analyzer 224 recalculates the average values 
for a and p, as well as the mean deviation values for a and p in step 342. In step 344, the statistical analyzer 224 

40 determines whether the deviation of a is less than 0.05. If the deviation ol a is less than 0.05, the result of decision 
344 is YES, and the system ends the process in step 352. If the deviation ol a is not less than 0.05, the result of decision 
344 is NO. In that event, the statistical analyzer 224 moves to decision 346 where it tests whether the deviation of p 
is less than 0.05 and the deviation of a is greater than 0.1. If the deviation of p is less than 0.05 and the deviation of 
a is greater than 0.1, the result ol decision 346 is YES. In that event, in step 350, the statistical analyzer 224 sets the 

45 average value ol a equal to the average value for p and ends the measurement process in step 352. If the result of 
decision 346 is NO. in step 348, the statistical analyzer 224 eliminates the farthest outlying values of a and p. Following 
step 348, the statistical analyzer 224 loops back to step 326 where it recalculates average values lor a and p and mean 
deviation values for a and p. This process is repeated until the deviation of a is less than 0.05. Following this process, 
the value of a and p may be displayed. The calculated value for a is also used in a conventional manner to determine 

so the S p 0 2 . This process is continually repeated for each new determination of an (a, P) pair. Thus, the value lor a is 
determined based on the statistical analysis of a plurality of previous (a, P) pairs. 

In one previously described embodiment of the system 200 the first and second data windows 302 and 304 (see 
Figure 11B) were derived using the AC component of the digitized signals 184 and 186 lor the first data window and 
the derivative of the digitized signals 184 and 186 lor the second data window. The example above utilized Equation 

55 21. which involves a number of integrals. However, as those skilled in the art will recognize, these calculations may 
beperformed on the individual data points by summing various values. This is illustrated in the flow chart of Figures 
1 3A and 1 3B. At the start 360. no data measurements have been made. The system 200 initializes the a and p arrays, 
which are stored in the a value storage area 21 8 (see Figure 9) and the p value storage area 220, respectively. In step 
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364, the system 200 retrieves 41 7 data points from the first and second data storage areas 204 and 206, respectively. 
The system 200 then enters the measurement loop. In step 368, the system gets an additional 100 data points. It 
should be noted that measurements are always done on the most recent set of 417 data points. In step 370, the system 
200 shifts the data to use the most recent 417 data points. In 6tep 372, the system 200 filters the data in the first and 

* second storage areas 204 and 206 using a 1 6 point average. This type of data averaging is well known to those skilled 
in the art and will not be described herein. 

In step 374, the mathematical processor 210 (see Figure 9) calculates a number of different values using the 417 
data points in the first and second data storage areas 204 and 206. These calculations include the mean, AC compo- 
nent, and derivative of each of the digitized signals 1 84 and 1 86. In addition, the mathematical processor 2 1 0 calculates 

io the sums of the square of the AC components, as previously described, the sum of squares of the derivative of the 
digitized signals 184 and 186, as well as the sums of the product of the AC components of the digitized signals and 
the sum of the product of the derivative of the digitized signals. The technique used to determine the mean value, the 
AC component, and the derivative has been previously described. In step 376, the window generator 214 generates 
the two data windows, as previously described. The first data window 302 (see Figure 11 B) is generated using the AC 

ts component of the digitized signals 1 84 and 1 86, while the second data window 304 is generated using the derivative 
of the digitized signals. < 

In step 378, the calculator 216 (see Figure 9) calculates values for a and p using the following equation: 

(a.M=t-*<»: 4ae),a y (23) 



where 
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a = ,^..Re dAC ..R AC . £ .R AC 2 .I^p.f).b = E 

XRed AC .|R AC . 



It should be noted that equation (23) is equivalent to equation (21 ) where the first data window 302 (see Figure 11 B) 
35 is the derivative of the digitized signals 184 and 186, while the second data window 304 is the AC component of the 
digitized signals. 

As previously discussed, it is not initially known which of the values corresponds to a and which value corresponds 
to p. However, in step 382, the system changes the order such that a is less than p to correspond with Equation (22). 
In decision 384, the statistical analyzer 224 (see Figure 9) determines whether the values for a and p are reasonable 

40 from a statistical point of view Specifically, in decision 384, the statistical analyzer 224 determines whether the max- 
imum value for the digitized signal 184 minus the minimum value for digitized signal is greater than one-third of the 
mean value for the digitized signal 1 84. In addilion, the statistical analyzer 224 determines whether the maximum value 
for the digitized signal 1 B6 minus the minimum value for that digitized signal is greater than one-third of the mean value 
for the digitized signal 186. If either of the maximum minus minimum values exceeds one-third of the mean value for 

45 the digitized signals 184 and 186, respectively, the result of decision 384 is YES, and the system 200 discards the 
current measurement in step 385, and uses the previous value for a in place of the current value. 

If the result of decision 384 is NO, or upon completion of step 385, the system 200, in step 386. the statistical 
analyzer scales the value of a and p by multiplying the measured values by the mean values. As previously discussed, 
this effectively removes the DC component of the digitized signals 184 and 186. In step 388, the system 200 shifts the 

so a and p arrays in the a value storage area 218 (see Figure 9) and the p value storage area 220 to contain the 1 2 most 
recent values. In step 390, the system 200 uses the analysis described in Figures 12A and 12B to calculate displayed 
values lor a and p. In step 392, the system calculates the oxygen saturation from the value of a in a conventional 
manner. In step 394, the system 200 displays the S p 0 2 value. Following the display of the S p O s value in step 394, the 
system 200 returns to step 368 and retrieves the next 100 data points. In this embodiment, the system 200 performs 

55 a new measurement approximately once per second. This illustrates a technique used to determine the value for a 
using two data windows. 

Alternatively, the system 200 can generate a plurality of windows and use other statistical techniques to derive a 
more accurate estimate for the value of a. A set of a values is generated for each of the plurality of data windows and 
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the final a value is selected from this set. The statistical analyzer 224 applies conventional statistical analysis techniques 
to the plurality of windows to derive the final a value. For example, the system 200 can use the first, second and third 
data windows 306, 308, and 310 (see Figure 11 C) to generate three different (a, P) pairs. 

Under ideal conditions, an intersection (not 6hown) of all three data windows 306, 308. and 310 on a graph would 
s provide the correct values for a and p. However, under less than ideal conditions, it is unlikely that there will be a single 
intersection between the curves generated by the data windows 306, 308, and 310. The statistical analyzer 224 will 
determine the most accurate a value in a manner described below. If the principles of the system 200 are extended to 
a large number of data windows, a more accurate statistical analysis may be performed to determine the values of a 
and p. For example, if the data sample 300 (see Figure 11 A) is divided into 100 data windows (not shown) the system 
io 200 would generate 100 (a, p) pairs. 

The calculator 216 (see Figure 9) uses equation (21 ) to solve for a and p in each of the 1 00 data windows. Thus, 
the calculator 216 determines 100 different values for the (a. P) pairs. The results of each of the calculations is stored 
in the a value storage area 218 as data bins. The use of bins is well known in the art and will only be described briefly 
herein. 

is The physiological range for a is subdivided into a number of ranges, designated as bins. Each time the calculator 
216 determines a value for a, the corresponding bin is incremented. Thus, the a value storage area 218 contains a 
frequency distribution of a values. This is illustrated in Figure 14 wherein the physiological range of values for a is 
subdivided into a plurality of bins 400. Each time the calculator 216 determines a value for a, the bin 400 corresponding 
to that value of a is incremented. This is indicated by an - X" in the bins 400. The bin 400 with the greatest number of 

20 occurrences can be selected as the best estimate for the value of ^Alternatively, the statistical analyzer 224 (see 
Figure 9) can analyze the a bin data and determine the value for a based on the results of the statistical analysis. For 
example, the statistical analyzer 224 can calculate the standard deviation of the frequency distribution of data in the 
a value storage area 218 and select the value for a based on the standard deviation (S.D.). In one embodiment, the 
statistical analyzer 224 selects an average value for a based on ±1 S.D. Other forms of statistical analysis well known 

25 to those of ordinary skill in the art may also be used. The present invention is not limited by the specific type of statistical 
analysis performed on the data. It should also be noted that the bin width illustrated in Figure 14 is provided as an 
example only. The bin width may be selected to provide the desired accuracy in determining the value for a. For 
example, the bin width illustrated in Figure 1 3 is 0. 1 However, if the physiological range for the value of ocean be divided 
into bin widths of 0.01 to provide a more accurate estimate for a. The bin width and the nature of the statistical analysis 

30 can be altered to suit the particular application in addition, the bin width and statistical analysis can be altered depending 
on the computing power ol the system 200. 

The method used by the system 200 for multiple data window pairs is illustrated in the flow chart of Figure 14. At 
the start 410, the sensor 2 (see Figure 4) has been placed on the subject's finger tip. In step 412, data is collected 
from the subject in the manner previously described. It should be noted that the ADC 124 provides digitized samples 

35 for each of the measured signals. R(t) and r(t). The digitized signals 184 and 186 (see Figure 9) are stored in the first 
data storage area 204 and the second data storage area 206, respectively. It should be noted, that while Figure 11 A 
illustrates a single data sample 300, the system 200 generates a data sample 300 for each of the digitized signals 184 
and 186. In step 414, the mathematical processor 210 (see Figure 9) calculates the derivatives, averages and the like, 
as described above. I n step 41 8. the window generator 21 4 generates the eight pairs of data windows described above. 

40 The calculator 216 solves for the a and p values in step 422 for all combinations of data windows or for selected 
combinations ol data windows. In step 424 the a values and p values determined by the calculator 216 are stored in 
the a value storage area 21 8 and the p value storage area 220, respectively. In step 426, the system 200 sets a equal 
lothe bin having the greatest number of occurrences. As discussed above, the statistical analyzer 224 may also perform 
statistical analyses on the data in the a value storage area 218. In step 428, the system 200 calculates the SpC 2 value 

45 in a conventional manner using the look-up table 14 (see Figure 4). The system 200 ends the measurement process 
at 430. 

As discussed above, the calculator 216. in step 422, calculates the a and p values for each of the generated 
windows. However, as discussed above, a large number of data windows can be generated. For each member of a 
pair of windows, the calculator 2 1 6 can determined an (a, P) pair using either the AC components, the measured signals 

so themselves, or the derivatives ol the measured signals. The calculator 216 can calculate a and p values for each 
possible permutation of windows. For example, the calculator 216 can calculate an (a, p) pair using the derivative 
values in window I and the AC values in window 1 , or the derivative values in window I with the AC values in window 
2. Thus, for each of the eight pairs of windows, the calculator 216 can generate 6 (a, p) pairs (see Table 1), resulting 
in 48 different (a, P) pairs. The 48 (a, P) pairs will generally be different. The binning process described above is one 

55 technique used to select the proper value lor a. However, those skilled in the art ol statistical analysis, will readily 
recognize that other techniques may be used to determine the best estimate lor the value of a. For example, sliding 
windows could be used instead of fixed bins. Alternatively, variable bin widths could be used. In yet another alternative, 
the system 200 subtracts the p count from the a count for each of the bins 320 (see Figure 12). Thus, the present 
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invention is not limited by the specific form of the statistical analysis. Once the value of a has been determined, the 
S p 0 2 is determined in a conventional manner. In addition, the peak S p O z value may be determined in a conventional 
manner using the S p 0 2 peak detector 136. 

It should be noted that the principles of the present invention may be extended beyond the measurement of blood 

5 oxygen saturation. For example, a third light source (not shown) may be added to produce a third wavelength in the 
sensor 2 (see Figure 4). Three ratios of light intensities (e.g., ratio of light source one to light source two, ratio of light 
source one to light source three, and ratio of light source two to light source three) can be derived from the light detector 
10. The three ratios can be independently used to derive both the arterial oxygen saturation and the arterial carboxy- 
hemoglobin saturation period. The present invention is also not limited solely to the use of optical sensors. Electrical 

10 sensors may derive physiological signals that can be processed according to the principles of the present invention. 
For example, electrical sensors can be used to derive a noise-free ECG signal. The electrical sensors each derive an 
ECG signal and the ratios of the ECG signals may be used to derive a noise-free version of the true ECG signal. 

In addition, the present invention may be used with a single sensor to measure a physiological parameter such as 
blood pressure. If the physiological parameter of interest is constant over the course of the measurement period, the 

is data sample may be divided into a plurality of data windows and the principles of the present invention applied to pairs 
of data windows. 

In operation, many of the components described above may be incorporated into a digital signal processor and/ 
or a digital computer. The programming details of the digital signal processor and computer are well known to those 
of ordinary skill in the art and need not be discussed herein. 
20 it is to be understood that even though various embodiments and~advantages of the present invention have been 

set forth in the foregoing description, the above disclosure is illustrative only, and changes may be made in detail, yet 
remain within the broad principles of the invention. Therefore, the present invention is to be limited only by the appended 
claims. 

25 

Claims 

1. A system for the enhancement of physiological signals for the measurement of blood oxygen in a subject, the 
system comprising: 

30 

first and second light sources to direct light toward the subject, said first and second light sources producing 
light of first and second wavelengths, respectively; 

a light detector positioned to detect first and second light signals alter interacting with the subject and to 
generate signals indicative of an intensity of said first and second detected light signals, said first detected 
35 signal having a first portion arising from light transmitted from said first source and a second portion arising 

from a first interference light source, said second detected signal having a first portion arising from light trans- 
mitted from said second source and a second portion arising from a second interference light source; 
a window generator to generate first and second data windows derived from data comprising said first and 
second detected signals; 

40 a storage location containing a mathematical relationship of said first and second portions of said first and 

second detected signals and a ratio of said first portion of said first detected signal to said first portion of said 
second detected signal; and 

a calculator using said mathematical relationship to generate first and second functions using said first and 
second data windows, respectively, said first and second functions having a point of intersection indicative of 
45 said ratio. 

2. The system of claim 1 wherein said ratio is indicative of blood oxygen saturation in the subject, the system further 
including a look-up table containing data relating said ratio to said blood oxygen saturation. 

so . 3, The system of claim 1 wherein said first data window is a derivative of said firsl and second detected signals. 

4. The system of claim 1 wherein said first data window is an alternating current (AC) component of said first and 
second detected signals. 

ss 5. The system of claim 1 wherein said first data window is derived from a group of signals comprising said first 
detected signal, a mathematical derivative of said first detected signal, said second detected signal, a mathematical 
derivative of said second detected signal, a combination of said first and second detected signals, and a mathe- 
matical derivative of said combination of said first and second detected signals. 
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6. Ttie system of claim 5 wherein said window generator calculates an average value for each of said group of signals 
and data less than said average for a particular one of said group of signals is placed in said first data window and 
data greater than or equal to said average for said particular one of said group of signals is placed in said second 
data window. 

5 

7, The system of claim 5 wherein said window generator calculates an average value for each of said group of signals 
and data less than said average for a first one of said group of signals is placed in said first data window and data 
greater than or equal to said average for a second one of said group of signals is placed in said second data window. 

io 8. The system of claim 1 wherein said mathematical relationship has the following form: 

jRf(t)-p jRi(t)ij(t) 



'* ai = Ft (P) = 1 



so 



jRi(0ri(t)-pfr?(t) 



where a is said ratio, i = 1 to 2 for said first and second data windows, respectively, Rj(t) is said first detected signal 
in the ith data window, r(t) is said second detected signal in the tth data window, and p represents a ratio variable 
of a portion of 6aid first detected signal to a corresponding portion of 6aid second detected signal. 

25 9. The system of claim 8 wherein F 1 (p)=F 2 (P), said calculator solving F( for ratio a 

10. A system for the enhancement of physiological signals for the measurement of blood oxygen in a subject, the 
system comprising: 

30 fi rs t and second light sources to direct light toward the subject, said first and second light sources producing 

light of first and second wavelengths, respectively; 

a light detector positioned to detect first and second light signals after passage through the subject and to 
generate signals indicative of an intensity of said first and second detected light signals, said first detected 
signal having a first portion arising from tight transmitted from said first source and a second portion arising 
35 from a first interference light source, said second detected signal having a first portion arising from light trans- 

mitted from said second source and a second portion arising from a second interference light source; 
a window generator to generate a plurality of data windows derived from data comprising said first and second 
detected signals; 

a storage location containing a mathematical relationship of said first and second portions of said first and 
40 second detected signals and a ratio of said first portion of said first detected signal to said first portion of said 

second detected signal; 

a calculator using said mathematical relationship to generate a plurality of functions from said plurality of data 
windows, respectively, each of said plurality of functions having a solution value indicative of said ratio; and 
a statistical analyzer to determine a final value for a ratio of said first portion of said first detected signal to 
45 said first portion of said second detected signal based on said plurality of solution values. 

1 1 . The system of claim 1 0 wherein said ratio is indicative of blood oxygen saturation in the subject, the system further 
including a look-up table containing data relating said ratio to said blood oxygen saturation. 

so 1 2. The system of claim 1 0 wherein a first of said plurality of data windows is derived from a group of signals comprising 
said first detected signal, a mathematical derivative of said first detected signal, said second detected signal, a 
mathematical derivative of said second detected signal, a combination of said first and second detected signals, 
and a mathematical derivative of said combination of said first and second detected signals. 

55 13. The system of claim 12 wherein said window generator calculates an average value for each of said group of 
signals and data less than said average for a particular one of said group of signals is placed in said first data 
window and data greater than or equal to said average for said particular one of said group of signals is placed in 
a second of said plurality of data windows. 
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14. The system of claim 12 wherein said window generator calculates an average value for each of said group of 
signals and data less than said average for a first one of said group of signals is placed in said first data window 
and data greater than or equal to said average for a second one of said group of signals is placed in a second of 
said plurality of data windows. 

1 5. A system for the enhancement of physiological signals representative of a physiological phenomenon in a subject, 
the system comprising: 

a sensor positioned in proximity with the subject to detect physiological signals and to generate signals indic- 
ative of said detected physiological signals, each of said detected signals having a first portion arising from 
the physiological phenomenon and a second portion arising from an interference source; 
a window generator to generate first and second data windows derived from said detected signals; 
a storage location containing a mathematical relationship of said first and second portions of said detected 
signals; and 

a calculator using said mathematical relationship to generate first and second functions using said first and 
second data windows, respectively, said first and second functions having a solution value indicative of the 
physiological phenomenon. 

1S. The system of claim 15 wherein said solution value is a point of intersection in said first and second functions. 

1 7. The system of claim 1 5 wherein said solution value is statistical measure based on said first and second functions. 

18. The system of claim 15 wherein said solution value is a ratio indicative of blood oxygen saturation in the subject, 
the system further including a look-up table containing data relating said ratio to said blood oxygen saturation. 

1 9. The system of claim 1 5 wherein said first data window is a derivative of said first detected signal. 

20. The system of claim 1 5 wherein said first data window is an alternating current (AC) component of said first detected 
signal. 

21. The system of claim 15 wherein the physiological signals are indicative of blood carbon monoxide levels, the 
system further comprising first, second, and third light sources generating first, second and third different wave- 
lengths of light, and a detector to detect light from said first, second, and third light sources alter interaction with 
the subject, said window generator generating first, second, and third data windows derived from said first, second 
and third detected signals, said mathematical relationship also including a first ratio of said first portion of said first 
detected signal to said first portion of said second detected signal and second ratio of said first portion of said first 
detected signal to said first portion of said third detected signal, said calculator using said mathematical relationship 
to generate first and second functions using said first and second data windows, respectively, said first and second 
functions having a solution value indicative of said first ratio related to blood oxygen saturation, said calculator 
using said mathematical relationship to generate a third function using said first and third data windows, respec- 
tively, said second and third functions having a solution value indicative of said second ratio related to blood car- 
boxyhemoglobin saturation. 

22. A method using a computer for the enhancement of physiological signals tor the measurement of blood oxygen 
in a subject, the method comprising the steps of: 

directing first and second light sources toward the subject, said first and second light sources producing light 
of first and second wavelengths, respectively; 

detecting first and second light signals after passage through the subject and generating signals indicative of 
an intensity of said first and second detected light signals, said first detected signal having a first portion arising 
from light transmitted from said first source and a second portion arising from a first interference light source, 
said second detected signal having a first portion arising from light transmitted from said second source and 
a second portion arising from a second interference light source; 

generating a plurality of data windows derived from data comprising said first and second detected signals; and 
using a mathematical relationship ol said first and second portions of said first and second detected signals 
and a ratio of said first portion of said first detected signal to said first portion o( said second detected signal 
to generate a plurality of functions from said plurality of data windows, respectively, said plurality of functions 
having a solution value indicative of said ratio. 
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23. The method of claim 22 wherein said plurality of functions comprises first and second functions, said first and 
second functions having an intersection indicative of said ratio. 

24. The method of claim 22 wherein said plurality of functions comprises first, second, and third functions, said ratio 
5 being statistically derived based on said solution values for said first, second, and third functions. 

25. The method of claim 22 wherein said ratio is indicative of blood oxygen saturation in the subject, the method further 
including the step of using a look-up table containing data relating said ratio to said blood oxygen saturation to 
determine the blood oxygen level in the subject. 

10 

26. The method of claim 22 wherein a first of said plurality of data windows is a derivative of said first of said plurality 
of detected signals. 

27. The method of claim 22 wherein a first of said plurality of data windows is an alternating current (AC) component 
is of said first of said plurality of detected signals. 

28. The method of claim 22 wherein af irst of said plurality of data windows is derived from a group of signals comprising 
said first detected signal, a mathematical derivative of said first detected signal, said second detected signal, a 
mathematical derivative of said second detected signal, a combination of said first and second detected signals, 

20 and a mathematical derivative of said combination of said first and second detected signals. 

29. The method of claim 22. further including the steps of calculating an average value for each of said group of signals 
and placing data whose value is less than said average for a particular one of said group of signals in a first of 
said plurality of data windows and placing data whose value is greater than or equal to said average for said 

25 particular one of said group of signals in a second of said plurality of data windows. 

30. The method of claim 22, further including the steps of calculating an average value for each of said group of signals 
and placing data whose value is less than said average for a particular one of said group of signals in a first of 
said plurality of data windows and placing data whose value is greater than or equal to said average for a second 

30 one of said group of signals in a second of said plurality of data windows. 

31. A method using a computer for the enhancement of physiological signals in a subject, the method comprising the 
steps of; 

35 directing a plurality of light sources toward the subject, each of said plurality of light sources producing light 

of a different wavelength; 

detecting a plurality of light signals after passage through the subject and to generate signals indicative of an 
intensity of said plurality of detected light signals, each of said plurality of detected signals having a first portion 
arising from light transmitted from a corresponding one of said plurality of light sources and a second portion 

40 arising from a plurality of interference light sources: 

generating first and second data windows derived from said plurality of detected signals; and 
a storage location containing a mathematical relationship of said first and second portions of said plurality of 
detected signals and a first ratio of said first portion of a first of said plurality of detected signals to said first 
portion of a second of said plurality of said detected signals: and 

45 using said mathematical relationship to generate first and second functions using said first and second data 

windows, respectively, said first and second functions having a solution value indicative of said first ratio. 

32. The method of claim 31 wherein said solution value is an intersection of said first and second functions. 
so 33. The method of claim 31 wherein said first ratio is statistically derived based on l said solution value. 

34. The method of claim 31 wherein said first ratio is indicative of blood oxygen saturation in the subject, the method 
further including the step of using a look-up table containing data relating said first ratio to said blood oxygen 
saturation lo determine the blood oxygen level in the subject. 

55 

35. The method of claim 31 wherein said first data window is a derivative of said first detected signal. 

36. The method of claim 31 wherein said first data window is an alternating current (AC) component of said first detected 
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signal. 

37. The method of claim 31 wherein said first data window is derived from a group of signals comprising said first 
detected signal, a mathematical derivative of said first detected signal, said second detected signal, a mathematical 
derivative of said second detected signal a combination of said first and second detected signals, and a mathe- 
matical derivative of said combination of said first and second detected signals. 

38. The method of claim 31 , further including the steps of calculating an average value lor each of said group of signals 
and placing data whose value is less than said average for a particular one of said group of signals in said first 
data window and placing data whose value is greater than or equal to said average for said particular one of said 
group of signals in said second data window. 

39. The method of claim 31 , further including the steps of calculating an average value for each of said group of signals 
and placing data whose value is less than said average for a particular one of said group of signals in said first 
data window and placing data whose value is greater than or equal to said average for a second one of said group 
of signals in said second data window. , 

40. The method of claim 31 wherein said first ratio is indicative of blood oxygen saturation in the subject, the method 
further including the step of using a look-up table containing data relating said ratio to said blood oxygen saturation. 

41 . The method of claim 31 wherein the physiological signals are indicative of blood carbon monoxide levels and said 
plurality of light sources comprises first, second, and third light sources, the method further including the steps of 
generating first, second, and third data windows derived from said first, second and third detected signals, respec- 
tively, said mathematical relationship also including a second ratio of said first portion of a first of said plurality of 
detected signals to said first portion of a third of said plurality of said detected signals, and using said mathematical 
relationship to generate first and second functions using said first and second data windows, respectively, said 
first and second functions having a point of inte rsection indicative of said first ratio related to blood oxygen saturation 
and using said mathematical relationship to generate a third function using said third data windows, respectively, 
said second and third functions having a point of intersection indicative of said second ratio related to blood car- 
boxyhemoglobin saturation. 

42. A method for the enhancement of physiological signals representative of a physiological event in a subject, the 
system comprising: 

positioning a plurality of sensors in proximity with the subject to detect a plurality of physiological signals and 
to generate first and second physiological signals indicative of the physiological event, said first and second 
physiological signals each having a first portion resulting form the physiological event and a second portion 
arising from interference: 

generating first and second data windows derived from said first and second generated signals; 
using a mathematical relationship of said first and second portions of said first and second generated signals 
to generate first and second (unctions using said first and second data windows, said first and second functions 
having a solution value indicative of the physiological event. 
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